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Abstract — We study the ergodic capacity of a fre- 
quency-selective Rayleigh fading channel with corre- 
lated scattering, which finds application in the area of 
UWB. Under an average power constraint, we consider 
a single-user, single-antenna transmission. Coherent 
reception is assumed with full CSI at the receiver and 
no CSI at the transmitter. We distinguish between 
a continuous- and a discrete-time channel, modeled 
either as random process or random vector with generic 
covariance. As a practically relevant example, we exa- 
mine an exponentially attenuated Ornstein-Uhlenbeck 
process in detail. Finally, we give numerical results, 
discuss the relation between the continuous- and the 
discrete-time channel model and show the significant 
impact of correlated scattering. 

I. Introduction 

Due to the increasing role of wireless communication 
it is important to determine the maximum achievable in- 
formation rates over multipath fading channels. Assuming 
an ergodic fading process and sufficiently relaxed decoding 
constraints, such that fluctuations of the channel strength 
can be averaged out, then the ergodic capacity is a suitable 
performance measure. It basically represents the average 
over all instantaneous channel capacities [1], [2]. 

In this paper, we examine the ergodic capacity of a 
single-user, single-antenna channel with full channel state 
information (CSI) at the receiver. The assumption of 
coherent reception is reasonable if the fading is slow in 
the sense that the receiver is able to track the channel 
variations. The transmitter has no CSI but knows the 
statistical properties of the fading. Further, we imply 
an average power constraint on the channel input and 
Rayleigh-distributed fading. 

Under the above constraints, the ergodic capacity was 
investigated for the case of flat Rayleigh fading, e.g., in 
[3], [4], [5]. The recent interest in ultra- wideband (UWB) 
technologies makes it important to examine the capacity 
also of frequency-selective fading channels. As related 
information-theoretic work, consider for instance [6], [7], 
where in [6] a model with two scattered components was 
studied. In [7] a multi-antenna system was considered. 

This research was supported by the Deutsche Forschungsgemein- 
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which can be adapted to serve as model for a single- 
antenna, frequency-selective fading channel. However, in 
either case uncorrelated scattering is assumed, which does 
not necessarily apply to UWB channels, as repeatedly 
documented in the literature [8], [9], [10]. This is one 
basic difference between conventional and UWB channels. 

In this paper, we study the ergodic capacity of frequency- 
selective fading channels with correlated scattering. To the 
best of our knowledge this has not previously been exa- 
mined in the literature. We consider models appropriate 
to characterize small-scale fading effects, i.e., fading due to 
constructive and destructive interference of multiple signal 
paths. We assume the small-scale fading to be Rayleigh 
distributed, which is not standard in UWB channel mode- 
ling [9]. However, detailed statistical evaluation of mea- 
surements in [10] support this assumption also for UWB. 
Furthermore, in [11] it is shown that the capacity of the 
Rayleigh fading channel provides a tight approximation, 
even if other fading statistics are employed. 

We distinguish between a continuous- and a discrete- 
time channel, where the channel impulse response (CIR) 
is either modeled as random process or random vector with 
generic covariance. Since both models are important, they 
are treated in parallel. The former is more general and 
more suitable to analysis, whereas the latter is more ade- 
quate for computer simulations and parameter estimation 
from measured data. Note, with the discrete approach we 
model equidistant samples of the CIR rather than variable- 
distant physical paths as in [12]. Modeling the sampled 
impulse response better describes the effective channel and 
is considered more robust since only aggregate physical ef- 
fects need to be reflected [10]. This is particularly relevant 
where frequency-selective propagation phenomena occur. 

For the continuous-time Rayleigh fading channel, we 
examine a detailed example with special covariance. We 
utilize an exponentially attenuated Ornstein-Uhlenbeck 
process being mathematically tractable and capturing 
an exponential power decay, which is common in UWB 
channel modeling [9], [12]. Additionally, it incorporates 
exponentially correlated scattering as measured in [8] . 

This paper is organized as follows: Section II specifies 
the channel models, and in Section III we derive respective 
expressions for the ergodic capacity. In Section IV we ana- 
lyze the example of an exponentially attenuated Ornstein- 



Uhlenbeck process, give numerical results, and discuss the 
relations between the continuous- and the discrete-time 
model. Section V finally concludes the paper. 

The following notation is used. The operator E[-] de- 
notes expectation, j is the imaginary unit, and X is the 
complex conjugate of X. We use the abbreviations R- 
integral for Riemann integral and i.i.d. for independent 
and identically distributed. A wide-sense stationary ran- 
dom process is referred to as stationary process. We define 
the sets Z^:={0, ...,A'- 1}, Kefi, and W:=(-f , -f ) , 
W > 0. Further, X = (Xk), kG 'Zj^, represents a column 
vector of size K with components Xk. Matrix notation is 
equivalent using two indices. 

II. Channel models 

As general fading multipath channel model we consider 

a linear time-varying system with equivalent lowpass im- 
pulse response (i?T,t) being a complex random process in 
the time variable tG'R and the delay parameter rSR. 
Then a realization h{T,t) := Hr^ti^) is the channel re- 
sponse at time t due to an impulse at time t — t [1], [13]. 
Next, for fixed r we assume the channel to be invariant 
within coherence intervals of fixed length. Hence, we may 
consider the random process (ffx,n)) t G M, in the discrete 
time variable neZ. Further, we imply (i?T,ra) to be sta- 
tionary and independent for fixed t, which corresponds to 
the block Jading model. As a consequence of this major 
simplification we are able to drop the time index n and 
model the channel as random process (-ffr) in the delay 
variable r G R. Another widely- used assumption of uncor- 
related scattering, i.e., E[HtHti] = for r, r'elR with 
T ^ t' , does not necessarily hold [8], [9], [10]. Therefore, 
we assume correlated scattering, which is a substantial dif- 
ference to previous work. 

In addition to fading we assume additive white Gaussian 
noise (AWGN) at the receiver. Below, we distinguish bet- 
ween a continuous- and a discrete-time channel. Thus the 
noise is either modeled as complex white Gaussian process 
{Zt),t € R, with i.i.d. real and imaginary part, each of zero 
mean and power spectral density ^ or as complex white 
Gaussian process with i.i.d. real and imagi- 

nary part, each of zero mean and variance 

Next, we will specify stochastic properties of (Ht) to 
obtain a Rayleigh fading channel model for continuous and 
discrete time. 

A. Continuous- Time Rayleigh Fading Channel Model 

Let {Xt), (^r)j T € R, be real stationary i.i.d. Gaussian 
processes with zero mean and continuous covariance func- 
tion R. Let g be an R-integrable function, i.e., g{T) dr 
exists as improper R-integral. We define 

ff(T):=«(T)l[o^^)(T),reR, (1) 

with 1^ the indicator function being 1 ii t € A and oth- 
erwise. Then {Xr) := (X^<?(r)), {Yr) := {Yrgir)), tGR, 



are real i.i.d. Gaussian processes with zero mean and co- 
variance function 

R{t, t') = R{t - r')g{r)g{T'), r, r' G R. (2) 

The indicator 1^ is introduced since only delays t>0 
are meaningful and u is some suitable function modeling 
power decay over r. Thus, (Xr), (Yr) are attenuated, 
non-stationary versions of (Xr), (Yr) with Xr = Yr = 
for T < 0. Finally, the continuous-time Rayleigh fading 
channel model is defined as (Hr) := (Xr + jYr), t G R. 

Note, by now nothing is said about band- or time- 
limitation of the channel. Combining a stationary process 
with a decaying function allows us to take advantage of 
well-investigated random processes as described, e.g., in 
[14, Ch. 3.5/3.7] while capturing the decaying nature of 
measured CIRs. For normalization we will use the constant 
c:=2 R{t,t) dr, which represents the mean energy 
contained in (Hr). It is finite if g^ is R-integrable. This 
is obviously a reasonable assumption from a practical 
viewpoint. Further conditions of R being continuous and 
g being R-integrable we will motivate later. 

Uncorrelated scattering can also be included in the con- 
tinuous model utilizing a covariance function R in form of 
a Dirac delta distribution, i.e., R{f) — cS{t) for some c > 0. 
However, to rigorously derive the ergodic capacity in this 
case we need to extend the mathematical tools applied in 
this work, as will be briefly discussed in Section III. 

B. Discrete- Time Rayleigh Fading Channel Model 

For the discrete-time model we assume the channel to be 
band-limited to W. Then (Hr) can be sampled at delays 
r = -j^, / G Z, to obtain the complex random process (Hi) 
in the discrete delay variable Z G Z. Note, we have infinite 
expansion in delay due to band-limitation and Hi ~0 for 
/ < due to i?T = for T < 0. In the following we will refer 
to Hi as the l-th channel tap. Next, we approximate (Hi), 
Z G by a random vector H := {Hq, . . . , H^^i) of size 
L := [WTfil -\- 1. Here, L models the number of significant 
channel taps with Td being the channel delay spread. This 
practically feasible approximation is well-founded since 
CIRs are sufficiently close to zero for delays t > T^, which 
is mathematically captured by u in (1). Note, for flat 
fading we have L = I, for frequency-selective fading we 
have L>1, and for UWB we clearly have 1. Finally, 
we denote the i-dimensional complex random vector H 
as H — X -hjY, where X, Y are real i.i.d. Gaussian vectors 
with zero mean and covariance matrix 

r := (7»fc), Ilk ■= QikCTiCTk, i.keZ^. (3) 

Therein, 2af is the mean power of the l-th channel tap 
Hi with af:=E[Xf]=E[Yi^]=E[\Hi\^]/2. If we deflne 
Pre ■=Pim ■= (crf), ^G Z^, then p:=Pre +Pim IS thc mean 
power delay profile, which is related to u in (1). The coef- 
ficients Qik represent the normalized correlation between 
tap Hi and Hk- Uncorrelated scattering is included as a 
special case, where the covariance satisfies r = diag(j3/2). 



For normalization we set the mean power of H to 1, i.e., 
= ^ • Subsequently, we refer to the above defined 
model as the discrete-time Rayleigh fading channel. 



III. Calculation of Ergodic Capacity 

We now calculate the ergodic capacity for the defined 
continuous- and discrete-time Rayleigh fading channel un- 
der the general conditions specified in the beginning of 
the introduction. Unless stated otherwise, capacity expres- 
sions are given in [bits/s] for the continuous model and in 
[bits/s/Hz] for the discrete model. 

A. Capacity Formulae 

In the continuous case the ergodic capacity within the 
frequency band W is calculated by [6], [1] 

/ log2(l + a|^/P)d/ , (4) 

,J-W/2 



C = E 



where (Hf) := {J^^e-^^^'f'^ Hrdr), f eM., is the Fourier 
transform of the process {Ht) and a: 



NoW 



defines the 



average signal-to-noisc ratio (SNR). The type of integral 
to calculate C and Hj is a stochastic R-integral as defined, 
e.g., in [15, Ch. 3.4]. Clearly, its existence is necessary for 
these quantities to make sense. Expression (4) is valid 
assuming an information carrying complex envelope input 
signal band-limited to W with constant power spectral 
density and mean power constraint P. 

For the discrete-time channel model (4) is also appli- 
cable if the discrete-time Fourier transform (DTFT), i.e., 
Hf := Ylfjo Hie-^^^'f^/^ J € W, is used to calculate the 
spectrum of the channel vector H. However, with the dis- 
crete approach we aim at a numerically easy-to-compute 
model, preferably discrete in the frequency domain as well. 
Therefore, we approximate the DTFT by an Af-point DFT, 
i.e., evaluating the spectrum at N points. Thus we calcu- 
late Hn ■= H f\f^nW/Nj nG'Zj^, and obtain the complex 
random vector H := {Hq, . . . , Hn-i). This actually means 
we are dividing the spectrum into A'' flat, parallel sub- 
channels which corresponds to an OFDM-based system ap- 
proach with N sub-carriers [2, Ch. 5.3.3/5.4.7]. The er- 
godic capacity is then given by 



Cm = E 



JV-1 



-^log2(l + a|F„|2) 

n=0 



(5) 



where a again means average SNR, now a:=-^. Con- 
sidering the parallel channels in frequency, the AWGN 
process at the receiver translates into a complex zero 
mean Gaussian random vector with independent real and 
imaginary parts with N independent components each 
having variance The average power constraint of P 
on each discrete-time channel input symbol converts to 
NP on the set of sub-channels (per OFDM symbol). Note 
that 2^ ^ C as ^ oo. 

In the following we evaluate (4) and (5). 



B. Continuous- Time Rayleigh Fading Channel 

Theorem 1: The ergodic capacity (4) of the continuous- 
time Rayleigh fading channel is given by 
rW/2 

C = RS) exp j Eii j df, (6) 

where 

R{t, t') cos (27r(r-T')/) dr dr' , / G R, (7) 

with R as defined in (2) and Ei„(z) := c^*H-"'dt. In 
particular, Eii(2;) =-Ei(-2;) with Ei the exponential inte- 
gral [16, 5.1]. If R in (2) satisfies R{f) = cd{f), we have 
uncorrelated scattering and obtain 



^=R^exp(^)Eii(^), 



(8) 



where c= 2c 9'^{t) dr is the only parameter of the CIR 
that has an influence on the capacity. Expression (8) is well 
known and was already derived in [6]. 

Proof: Here we just provide the main parts of the 
proof that allow to understand the underlying methods. A 
complete proof containing omitted details is given in [17]. 
Relevant properties of stochastic R-integrals and of mean 
square calculus can also be found in [18, Ch. 2.1-2/8.1-2], 
[15, Ch. 3.3/3.4/3.6], [19, Ch. 8-4.]. 

(i) Stochastic R-integral: Let (?7^), r e R, be a complex 
random process with i!^[|J7T-|^] < oo. Then it can be shown 
that the stochastic R-integral Jj/ := 11^- dr exists if and 
only if E[UrUr'] is R-integrable, i.e., 



/OO pC 
-OO ^ —( 



E[UrUr']dTdT' 



(9) 



exists as improper R-integral. Then we obtain 

/OO 
E[Ur] dr. (10) 
-OO 

If we have another complex process (V^), rGlR, with 
£?[|Vr|^] < OO for which ly '■= dr exists, then it can 

further be shown that 

/OO /.OO 
/ E[UrVr']dTdT'. (11) 
-OO ./ — OO 

(a) Existence of (Hf): We use (9) to prove that the 
Fourier transform {Hf) = {f^^e-i^''f''HrdT), /sR, of 
the CIR (Hr), r e R, exists. Thus we set Ur := c-^^tt/t^^ 
and obtain (9) with E[UrUr>]=2R{T,T')e-^^''f^e^^''f''' , 
where R is the covariance given in (2). If R in (2) is such 
that R is continuous and g is R-integrable then (9) exists 
for all / gR and hence (Hf). Note, for {Hf) to exist, we 
can alternatively require R to be R-integrable. 

(Hi) Distribution of (Hf): It can be shown that any 
linear transformation of a Gaussian process is a Gaussian 
process. Thus (Hf) is a complex Gaussian process com- 
posed of two real Gaussian processes {Xf) := (Re[Hf]) and 
{Yf) := {lm[Hf]). The process (Xf) has zero mean, which 



follows from (10) with [7^ := Re[e-^'2'^/^i?^] and from 
E[Hr] — 0- Wc calculate the covariance of (Xf) using (11) 
with Ur as before and Vr '■=Re[e~^'^'^-^ Ht] resulting in 
E[XfXf,] = /,°° /,°° i?(T, r') cos {2n{Tf - r' f)) dr dr' for 
/, /' G R and R as in (2). Correspondingly, we show with 
Re[-] replaced by Im[-] that the process (Yf) has zero mean 
and identical covariance. Again we use (11) with Ut '■= 
Re[e-J'2'^^^iJ^] and Vr :=Iin[e-:>'^''f'^ Hrj to determine 
the cross-covariance between (Xf) and {Yf). We obtain 
E[XfYr] = XT JT ^(^' (Mrf - r'f)) dr dr' for 

/,/'g]R and E[XfYf,] = -E[Xf,Yf]. Thus the processes 
{Xf), {Yf) are identically distributed but not independent. 

(iv) Distribution of\Hf\^: From (iii) it follows for any 
/gR that Xf,Yj are i.i.d. Gaussian random variables 
having zero mean and variance E[Xj\ =d'^{f) with d^{f) 
as in (7). Thus \Hf\'^ =Xj + Yj is an exponentially dis- 
tributed random variable with probability density function 
qf{z) = (2a2(/))-i exp[(2a2(/))-i£], z > 0. 

(v) Existence and calculation ofC: The criterion (9) for a 
stochastic R-integral to exist and the properties (10), (11) 
are also valid for finite integration boundaries. Thus if we 
set Uf :=log2(l + a |-ff/P) then lu := /^/j 'V exists if 
and only if E[UfUfi] is R-integrable over Due to the 
finite integration boundaries it is sufficient to show that 
E[UfUf'] is continuous. This is equivalent to E[Uj] being 
continuous. The continuity of E[Uj] can either be shown 
directly or equivalently by showing that the process {Uf) 
is continuous in mean square. Now we can apply (10) and 
obtain C = E[Iu] = ^i^A '^f- We calculate E[Uf] by 
evaluating the integral E[Uf]~ /Q°°log2(l + az)qf{z)dz 
with Qf as in (iv). Finally, this yields (6). 

(vi) Uncorrelated scattering: The covariance R in (2) 
becomes R{t,t') =cS{t — T')g{T)g{T') for R{t) = cS{f). 
Plugging R into (7) and using the fundamental property 
of the Dirac distribution of h{T)S{T — t') dr = h{T') 
yields {f) = c g"^ {t) dr , independent of /. Then 
(6) immediately implies (8). Note, this is only a formal 
derivation. The introduced mathematical tools are not 
applicable to random processes with Dirac-type covariance 
functions, which do not even exist in the classical sense. 
A rigorous mathematical treatment involves stochastic 
differential equations. ■ 

C. Discrete-Time Rayleigh Fading Channel 

Theorem, 2: The ergodic capacity (5) of the discrete- 
time Rayleigh fading channel is given by 

CN = 7rWT)Il {2^) (2^) (12) 

n=0 

where 

L-l L-1 

^l = J2Yl (27r(i - k)n/N), n G Z^^, (13) 

i=0 fe=0 

with 7ife as in (3). If (3) satisfies r= diag((TQ, . . . , (y\_i), we 
have uncorrelated channel taps (uncorrelated scattering) 



and obtain 

C;v = 45yexp(l)Eii(l)=:a3, (14) 

independent of L, N, and the mean power delay profile 
(ctq, . . . , cr£_i)- Expression (14) is identical to the capacity 
of the flat Rayleigh fading case as derived in [3], [4]. 

Proof: Theorem 2 is similarly proved as Theorem 1 but 
with less effort. We may write the iV-point DFT H of the 
channel vector H using matrix notation. Let ^ := {(pni) 
with (^„, :=e-^'2"'/-^^, neZ^, /gZ^. then H = <i>H. The 
existence of H is now evident. Further, transforming a 
complex Gaussian vector this way yields a complex Gaus- 
sian vector [18, 7.5-2] with X := Re[.ff], Y :=lm[H] being 
real Gaussian vectors. We calculate means and covariances 
using the linearity of E[] and the independence of X = 
Re[-ff] and Y = lm.[H]. We obtain that X, Y are identically 
distributed with zero mean and covariance E[XnXn'] = 
Y.i=o I]fc=o 7ife cos (27r(m - kn')/N) for n,n' eZ and 
"fik as in ^). The cross-covariance is given by i?[A„y„'] = 
J2i=o J2k=o ^li-k {^^{in — kn')/N) for n, n' G Zjy with 
E[X„Y„,] ^ -E[Xn'Ya]. It follows for all n G that X„, 
Yn are i.i.d. Gaussian random variables with zero mean 
and variance E[X^j]=d-'^^ with ct^ as in (13). Therefore, 
\Hn\'^ = X^ -|- 1^ is again an exponentially distributed 
random variable with parameter {2a^)~^. Now we ex- 
change summation and expectation in (5) and evaluate 
E[log2{l + a\Hn\'^)] to obtain (12). Finally if we have 
uncorrelated scattering, i.e., F in (3) is diagonal, then 

= Ef=V = 5 forn e Zjv, and (12) implies (14). ■ 

IV. Example: Exponentially attenuated 
Ornstein-Uhlenbeck process 

Now, we consider an example for the continuous-time 

Rayleigh fading channel with special covariance. We use 
an exponentially attenuated Ornstein-Uhlenbeck process 
capturing an exponential power decay. This is a common 
assumption in UWB channel modeling and was justified 
by measurements [9], [12]. In addition, it incorporates ex- 
ponentially correlated scattering. Based on UWB channel 
measurements, a similar correlation model was used in [8]. 
Finally, we discuss relations to the discrete-time model. 

A. Definition 

A stationary Ornstein-Uhlenbeck process is a real Gaus- 
sian processes with zero mean and covariance fimction 
R{t — t') := ^e""!"^""^ I for r, t'gR with parameters 
a,d>0 [14, Ch. 3.7.2/3.7.3]. Using the notation of Sec- 
tion II- A wc set {Xr), {Yr) to be normalized Ornstein- 
Uhlenbeck processes with d:=2a and specify the attenu- 
ation function u in (1) as u{t) ■.= ^/bce~^'^ for tGR with 
parameters 6, c > 0, where c is the constant defined in Sec- 
tion II-A. We obtain the attenuated Ornstein-Uhlenbeck 
processes [Xr) = {Xrg{T)), {Yr) = {Yrg{T)), relR, repre- 
senting independent real and imaginary part of {Ht), each 
with covariance function 

R{ry) = ce-"l---'l6e-^(-+^')l[o,^)(r)l[o,^)(r'), (15) 



for r, r' e R. 

B. Analytical Calculations 

The expressions given in this subsection are derived in 
the Appendix in condensed form. 
Calculating (7) using (15) we get 

c(a + b) 



(16) 



(a + 6)2 + (27r/)2 

by elementary integration. One representation of the closed 
form solution of (6) is given by 



r — ^ \ " _i_ \ ' 

~ kipy n+i 

n=0 k=0 



^LU(3,WyA)Lj,{(3M (17) 



where is the generalized Laguerre polynomial of order k 
[20, 8.970], /3i :=27rV(ac(a + 6)), and p2 := (a + 6)V(47r2). 
Upper and lower bounds for (17) are given by 



Ce =Wlog2 (l + 2e-^aa'^{W/2)) 



ln(2) 



arctan 



W/2 



(18) 



where we have an upper bound for 9 = Q and a lower bound 
for = 7 with 7 being the Euler constant. 

In case the effect of the channel outside the band W 
is negligible, i.e., W is sufficiently large (depending on 
the parameters a, 6, c, a), then (17) can be closely ap- 
proximated by integrating (6) with (16) over entire R 
yielding 



:exp(/3i/32)r(|,/3i/32), 



(19) 



ln(2)Vft 

where T{^i,z):= e~*t^~^dt is the incomplete gamma 
function [20, 8.35o''.2]. 

Note, C~ in (19) is just an approximate expression for 
proper parameter ranges but not the capacity for infinite 
bandwidth, i.e., not limv^^oo C' since a and thus /3i de- 
pend on W as well. 

C. Numerical Results and Relation to the Discrete- Time 
Rayleigh Fading Channel 

Here, we give numerical examples for the previously de- 
rived expressions and show the connection to the discrete- 
time channel. To do so, we first define two constants for 
the continuous-time channel. Let eG (0, 1) be the portion 
of the channel within the frequency band W in terms of 
mean energy, i.e., £ := ^ '^'^(f) 'V > where c is defined 

in Section 11- A. If we fix i, then we get W = ^ tan(|e) by 
elementary integration. Further let e G (0, 1) be the portion 
of the channel within the time interval [0, T,;] in terms of 
mean energy, i.e., e := | Jq"^ R{s, s) ds. If we fix e, then we 
get Td = ^ ln(l — e) by elementary integration. 

Consider the continuous-time channel within the fre- 
quency band W. To represent this channel by the discrete- 
time model we sample the covariance function R of (15) 
over the range [OjTd] x [0,Td] by (^-spacing to get the 
covariance matrix F of (3). We normalize F to have channel 



vectors with unit mean energy as stated in Section II-B. 
In addition, we choose e to be close to 1, i.e., truncation in 
time is negligible. Note, this discrete version of the channel 
resembles the behavior of the continuous channel in time 
domain. However, the spectrum is different due to aliasing 
since the Ornstein-Uhlenbeck process is not band-limited 
because of (16). This, in turn, is negligible if e is close to 1. 
To obtain equality in frequency domain, we have to sample 
a lowpass-filtered version of the continuous channel, which 
of course has different behavior in time. Since we aim at 
equivalence in time, we consider the former approach. 
Finally, for the discrete and the continuous model to be 
comparable we have to set c:=W/£. 

As a numerical example, we set a:=b:=^, and e := 
0.998 resulting in Td = 6.2146, when normalized to sec- 
onds. In Example 1 (Fig. 1), we set e:= 0.998 leading 
to VF= 101.32, when normalized to [Hz], c= 101.52, and 
L = 630. In Example 2 (Fig. 2), we set £:= 0.800 leading 
to = 0.9796, c= 1.2245, and L = 7. Capacity expres- 
sions are considered as function of a and are given in 
[bits/s/Hz], where normalization to the respective band- 
width W is performed if required. The values of C/W 
are obtained by numerically evaluating (6) and Cjv is 
computed with N = 6300. As reference curves the AWGN 
channel capacity Cawgn = log2(l + oi) is plotted. Another 
reference is the capacity for uncorrelated channel taps 
(14), which is an upper bound for the correlated case due 
to Jensen's inequality. 

We observe differences between C^ and C/W which are 
due to the applied approximations, i.e., truncation in time, 
aliasing, and discretization of spectrum. They are prima- 
rily minor, particularly in Example 2, but increase with a. 
The distance to Cus is considerable in Fig. 1 but small in 
Fig. 2. This depends on the concentration of the spectrum 
within the considered band, which in turn is controlled 
by the degree of correlation (parameter o) and the de- 
lay spread (parameter b). The tightness of the bounds is 
parameter-dependent. Especially the lower bound in Fig. 2 
is very tight for a > 15dB. Finally, C is well approximated 
by Cr; in Example 1 for a < 15dB but is inappropriate in 
Example 2, since e is not close enough to 1. 

V. Conclusion 

In this paper, we determined the ergodic capacity of 
a frequency-selective Rayleigh fading channel with cor- 
related scattering. We considered a continuous- and a 
discrete-time channel and examined a detailed exam- 
ple incorporating exponential power decay and exponen- 
tially correlated scattering. Analytical, approximate, and 
bounding expressions as well as numerical results were pre- 
sented and the relation between continuous- and discrete- 
time models were discussed. The results illustrate sig- 
nificant differences between the capacities for correlated 
and uncorrelated scattering. Future work includes for 
example: considering other stationary processes, assuming 
non-perfect CSI, further elaborating (17), analyzing the 
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Fig. 1. Ergodic capacity Example 1: a = b = ^,6 = 0.998, i = 0.998. 
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Fig. 2. Ergodic capacity Example 2:a = b=^,e = 0.998, i = 0.800. 

outage capacity, or estimating model parameters from 
measured data. 

Appendix 

A. Derivation of (17).- 

(i) Using (16) in (6) and Ei,{x)e- = -A_L„(x) 
[21, 5.11.1.4] we get C=j^ Y^Zo T^i^i ^it^ ■= 
I-w/2 LniPiih + P)) df by Lebesgue's theorem. 

(ii) Using the substitution s = /3i/^ and Ln{x + y) = 
YZ^^L-HxWXiy) [21, 4.4.2.3] then yields h = 
1^1 ^ Z]fc=o -^rt-fc(/5i/^2)^2 with the remaining integral 
h := Jo s 2Lf^^{s)ds. 

(iii) Finally, we evaluate the term I2 using s^L^(s) ds = 

[21, 1.14.3.4] to get (17). 

B. Derivation of (18); 

We simply use the inequahtics ln(l + ^) < Eii(a:;)e^ 
[22, eq. (13)] and Eii(a;)e^ < ln(l + ^) [16, 5.1.20] to re- 
place the integrand. Then elementary integration and the 
monotony of the integral yields (18). 

C. Derivation of (19).- 

(i) Using (16) and infinite integration boundaries in 
(6) we get = ^^^^ /~ ^ Eii(s + A/J^) ds by 



again substituting s = /3i/^. 

(ii) With/- ^Eii(^+.)d^'^ -';:X;r^V ^^^(.) 

[21, 2.5.3.14] and W_i^i{y) = y-^e'iT{\,y) [20, 
9.236.1, 8.359.3], where V„,a is the Whittaker's 
W-function [20, 9.220], we finally get (19). 
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